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Abstract. Josephson-junction systems have been studied as potential realizations of quantum 
bits. For their operation as qubits it is crucial to maintain quantum phase coherence for long 
times. Frequently relaxation and dephasing effects are described in the frame of the Bloch 
equations. Recent experiments demonstrate the importance of 1// noise, or operate at points 
where the linear coupling to noise sources is suppressed. This requires generalizations and 
extensions of known methods and results. In this tutorial we present the Hamiltonian for 
Josephson qubits in a dissipative environment and review the derivation of the Bloch equa- 
tions as well as systematic generalizations. We discuss 1/f noise, nonlinear coupling to the 
noise source, and effects of strong pulses on the dissipative dynamics. The examples illus- 
trate the renormalization of qubit parameters by the high-frequency noise spectrum as well as 
non-exponential decay governed by low-frequency modes. 



1. Introduction 

Quantum-state engineering requires coherent manipulations of suitable quan- 
tum systems. The needed quantum manipulations can be performed if we 
have sufficient control over the fields which couple to the quantum degrees 
of freedom, as well as the interactions. The effects of external noise sources 
have to be minimized in order to achieve long phase coherence times. In 
this tutorial we review the requirements for suitable physical realizations of 
qubits - with emphasis on Josephson circuits - and discuss ways to analyze 
the dissipative effects. 

As has been stressed by DiVincenzo (1997) any physical system, considered 
as a candidate for quantum computation, should satisfy the following five 
criteria: (i) One needs well-defined two-state quantum systems (qubits). (ii) 
One should be able to prepare the initial state of the qubits with sufficient 
accuracy, (iii) A long phase coherence time is needed, sufficient to allow for 
a large number (depending on details, say, > 10 4 ) of coherent manipulations, 
(iv) Control over the qubits' Hamiltonian is required to perform the necessary 
unitary transformations, (v) Finally, a quantum measurement is needed to 
read out the quantum information. 

The listed requirements may be satisfied by a system of spins - or quantum 
degrees of freedom which under certain conditions effectively reduce to two- 
state quantum systems - which are governed by a Hamiltonian of the form 

H = H ctv i(t) + H meas (t) + Hj. es . (1) 
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The first term, H ctr \(t), describes the control fields and interactions, 

ffdxi(t) = ~E E &)<<4- (2) 

i,a iy^j;a,b 

Here a a (with a,b = x, y, z) are Pauli matrices in the basis of states ||) = 
(j) and ||) = (°). Full control of the unitary quantum dynamics of individ- 
ual spins is achieved if the effective 'magnetic field' B l {t) can be switched 
arbitrarily at each site i. For most purposes control over two of the field 
components is sufficient, e.g., H cU \(t) = — | [B\ (t)a l z + B l x (t)a x ] . In 

i 

order to perform logic operations, e.g. for quantum computing, one also needs 
two-qubit operations. They can be controlled if the coupling, J^ b (t), between 
the qubits can be switched. Examples for a suitable coupling are an Ising 
zz-coupling or a spin-flip xy-coupling (see below). 

The measurement device and the residual interactions with the environment 
are accounted for by extra terms H meas (t) and H rcs , respectively. Ideally the 
measurement device should be switchable as well and be kept in the off- 
state during manipulations. The residual interaction H ICS leads to dephasing 
and relaxation processes. It has to be weak in order to allow for a series of 
coherent manipulations. 

A typical experiment involves preparation of an initial quantum state, switch- 
ing the fields B l (t) and coupling energies J^ b (t) to effect a specified unitary 
evolution of the wave function, and the measurement of the final state. 
The initial state can be prepared by keeping the system at low temperatures 
in strong enough fields, B\ 3> ksT, for sufficient time such that the residual 
interaction, H ICS , relaxes each qubit to its ground state, |f ). 
A single-bit operation on a selected qubit i can be performed, e.g., by turn- 
ing on the field B l x (t) for a time span r. As a result the spin state evolves 
according to the unitary transformation 

Ui (a) = exp — JL = I , 2 2 (3) 

where a = B 1 x t/H. By appropriate choice of the parameters an a = it- or an 
7r/2-rotation can be induced, producing a spin flip (NOT-operation) or (start- 
ing from the ground state ||)) an equal-weight superposition of spin states, 
respectively. Turning on B l z {t) for some time r produces another needed 
single-bit operation, a phase shift between ||) and ||), described by U\ {[3) = 
exp (iB l z T(j l z /2fi) where [5 = B 1 z t/H. With a sequence of x- and z-rotations 
any unitary transformation of the single-qubit state can be achieved. 
A two-bit operation on qubits i and j is induced by turning on the correspond- 
ing coupling J % a b (t). For instance, for the xy-coupling (a x a J x + cr l y cr 3 y ) the 
result is described, in the basis |T;fj), Itilj), l-Utj). by the unitary 
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Figure 1. a) A Josephson charge qubit in its simplest design formed by a superconduct- 
ing single-charge box. b) A qubit with the electromagnetic environment represented by the 
impedance Z(lo). 



operator 



ugh) 



( 1 \ 

cos 7 ism 7 
ism 7 cos 7 

Vo o o 1/ 



(4) 



with 7 = 2J 1 ^t/H. For 7 = tt/2 the operation leads to a swap of the states 
ITilj) and ||; Tj) ( an d multiplication by i), while for 7 = tt/4 it transforms 
the state into an entangled state (\Uij) + i 1-Utj))- 

Instead of the sudden switching of B l z x (t) and J l i(t), discussed above for 
illustration, one can use other techniques to implement single- or two-bit 
operations. For instance, ac resonance signals can induce Rabi oscillations 
between different states of a qubit or qubit pairs. Both switching and ac- 
techniques have been applied for Josephson qubits, e.g., by Nakamura et al. 
(1999) and Vion et al. (2002), Chiorescu et al. (2003), respectively. 
The coupling to the environment, described by H ICS , leads to dephasing and 
relaxation processes. In this tutorial we will first derive in Section 2 the 
proper form of H TCS for the case where a Josephson charge qubit is coupled 
to an electromagnetic environment characterized by an arbitrary impedance. 
In Section 3 we present a systematic perturbative approach, derive the Bloch 
equations and discuss their validity range. In particular, we find expressions 
for the relaxation and dephasing rates. In the following Section 4 we use the 
Bloch equations to study several problems of interest. Section 5 deals with 
extensions beyond the Bloch-Redfield description. 
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2. Dissipation in Josephson circuits 



2.1. Josephson qubits, the Hamiltonian and the dissipation 

Fig. 1 shows an example of a 'Josephson charge qubit' built from supercon- 
ducting tunnel junctions. All Josephson qubits have in common that they are 
nanometer-size electronic elements embedded into and manipulated by elec- 
trical circuits. The electromagnetic environment, which can not be avoided if 
we want to control the qubit, leads also to dissipation. This environment plays 
a crucial role in many contexts, e.g., the physics of the Coulomb blockade in 
tunnel junctions, for which the appropriate theory is reviewed by Ingold and 
Nazarov (1992). In the spirit of this analysis we will analyze the effects of the 
environment on Josephson qubits. 

We first consider a superconducting single-charge box shown in Fig. la in 
the absence of dissipation. We further choose parameters such that single- 
electron tunneling is suppressed and only Cooper-pair charges tunnel across 
the junction. This situation is realized if the superconducting gap A is larger 
than the charging energy scale Eq (defined below) and the temperature k^T. 
In this case superconducting charge box is described by the Hamiltonian 



Here n is the operator counting the number of excess Cooper pairs on the 
island relative to a neutral reference state. It is conjugate to 9, the phase 
difference of the superconducting order parameter across the junction. Charge 
quantization requires that physical quantities are 2ir -periodic functions of 6, 
for instance e t6 ', with e ld \n) = \n + 1). The 'gate charge' Q g = C g V g is 
controlled by the gate voltage, and Eq = e 2 /2(Cj + C g ), which depends 
on the junction and gate capacitances, is the single-electron charging energy 
scale. Interactions with charges in the substrate induced by stray capacitances, 
leads to additional contributions to the gate charge. 

Next we account for the dissipation due to electromagnetic fluctuations. They 
can be modeled by an effective impedance Z{w), placed in series with the 
voltage source (see Fig. lb) and producing a fluctuating voltage. The sum 
V g + SV(t) couples to the charge in the circuit. For an impedance with open 
leads the Johnson-Nyquist expression relates the voltage fluctuations between 
its terminals to Z(uS). In the present case the impedance is embedded in a cir- 
cuit, which further modifies the spectrum of voltage fluctuations. Neglecting 
for a moment the Josephson tunneling (i.e., treating the Josephson junction 
as a capacitor) we obtain from the fluctuation-dissipation theorem (FDT) 



SvW = 77 / dte lu}t ({5V{t),SV(0)}} = ReZ t (u) tua coth — — . (6) 




(5) 



fojj 




2k B T 
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Here Z t (uj) = [-tuC t + Z' 1 ^)] 1 and C t _1 = Cj l + C" 1 are the total 
impedance and capacitance as seen from the terminals of Z(ui), respectively. 
For many purposes it is sufficient to characterize the environment by its noise 
spectrum and the corresponding response function. In general one should 
take into account that the environment itself is also a quantum system with 
many degrees of freedom. As argued before (Caldeira and Leggett, 1983), in 
a generic setting one can think of a linear environment as of a large set of har- 
monic oscillators, each of which interacts weakly with the system of interest. 



2m a 



2 (l 



, and the qubit is coupled to the fiuctua- 
J2 a ^aX a - The Johnson-Nyquist relation (6) is reproduced 



I-e., -ffbath = J2a 

tion variable 5V = 

if one chooses the bath 'spectral density' as Jy (uj) = | J2a ml) ^{ UJ ~ LjJ a) = 
oj ReZ t (uj). The Hamiltonian of the whole system then reads 1 



H = H hox - 2en^ 5V + H hsith + -t- 



2en 



Cj 



(7) 



The last 'counter-term' is introduced to cancel the renormalization of the 
qubit's charging energy by the environment. The capacitance C t is defined by 
Jq°° dm Jy{uj)/uj = 7r/(2Ct). Because of the analytic properties of response 
functions it also equals C t = C t + lim^oo i/(uoZ(uj)). For a pure resistor 
Z(lo) = R it reduces to C t . The last three terms in (7) can be lumped as 



h — H hox + 



Pi m a ujl 
2m n 2 



Q \ 2 

. 2677- 



2.2. Reduction to a two-level system 



If we bias the superconducting charge box near a degeneracy point of the 
charging energy, e.g. at a voltage close to e/C t , and choose all other parame- 
ters (temperature, frequency and strength of control pulses) appropriate, only 
two charge states play a role, |n = 0) and |n = 1), which we denote as ||) 
and ||), respectively. The leakage to higher charge states can be kept at a 
low level (Fazio et al., 1999). In spin notation, the number operator becomes 
n = \(l — a z ), while cos 9 = \o x , and the effective spin Hamiltonian reads 

H = --B z {V g ) u z - -B x a x -^Xa z + H h&th , (8) 

1 Here we 'postulate' the model based on experience with dissipation in quantum systems. 
For the reader who is less familiar with these ideas we present in the appendix a more detailed 
derivation of the Hamiltonian for a particular model of the electromagnetic environment, 
namely a resistor modeled as a transmission line. 
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where X = 2e5V ,B X = E 3 , and 



c, 

B ' = c?Tc^- c ^ = iE °{ 1 - < T^ <9) 



The spectral density of the fluctuating field X is related to the power spectrum 
of voltage fluctuations by 

Sx{u) = \ ((Xl) + (X^)) = ( 2e g) 2 Sy(u) . (10) 

Here we defined {X*) = J dt e iu)t (X(t)X(0)}. For the calculation of many 
properties of the two-level system the knowledge of the symmetrized noise 
correlator, Sx, is sufficient. In general, however, the anti-symmetric part 
Ax{uj) = \ ([X (t) , X (0)]} ^ is needed as well. It is related to the response 
function of the bath x(t) = (i/h)9(t)([X(t),X(0)}), which describes the re- 
action of the bath to a perturbing force. Namely, if a perturbation H p = — fX 
is added to the Hamiltonian the bath responds as (X) w = x(a;)/(u;). The 
function Ax(u>) is related to the imaginary part of the response function, 
Ax{uj) = hx" (w), and in equilibrium the FDT fixes Sx(w) = Ax{uj) coth(frw/2kBT). 
To simplify the comparison with the literature on the Caldeira-Leggett model 
we include the pref actor in (10) into the definition of the spectral density, 
Jx(u) = (2eC t /Cj) 2 Jy(a;) = x"(w). The generic low-frequency behavior 
is a power-law, Jx(^) = 2nhaujQ~ s uj s , where ujq is a frequency scale and 
a the dimensionless strength of the dissipation. Of particular interest is the 
Ohmic case (s = 1), obtained if Z(u>) = R and, hence, ReZ t (uj) rs R, for 
u) <C (RCt)^ 1 . In this case we have Jx = 2nhauj and 

S x (oj) = 2TTh 2 aujcoth^-, (11) 

with 

and Rq = (2e) 2 /h is the (superconducting) resistance quantum. 



3. Dissipative dynamics 

3.1. B LOCH EQUATIONS 

The dissipative dynamics of spins has been the subject of extensive research 
in the context of the nuclear magnetic resonance (NMR). One of the main 
tools in this field are the Bloch equations. These kinetic equations were first 
formulated, on phenomenological grounds, by Bloch (1946) for the case of 
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nuclei in a magnetic field B = B\\ + 5Bj_(t), which is the sum of a strong 
static field By in the z-direction, and a weak, time-dependent, transverse per- 
turbation 5B± (t). The latter may be chosen to oscillate in resonance with the 
Larmor frequency B»/H to induce spin flips. The Bloch equations describe 
the dynamics of the magnetization in the 'anisotropic r-approximation' : 

±M = -BxM-±-(M z - M )z- ±-(M x x + M y y) . (13) 

The two relaxation times, T\ and T 2 , characterize the relaxation of the longi- 
tudinal component of the magnetization M z to Mq, which is the equilibrium 
magnetization in the static field B» , and the transverse components (M x , M y ) 
to zero, respectively. Eq. (13) describes the evolution of a two-level system in 
terms of its 'magnetization' M cc (a), and at the same time it describes the 
time evolution of the components of its density matrix, related to the magne- 
tization by (a z ) = poo — Pn and (fx + ioy) = Pio- Here and 1 denote the 
ground state, and excited state, |J, Z ), in the field B», respectively. 
Using the normalization, trp = poo + Pn = 1> one can rewrite Eq. (13) as 
equations of motion for the density matrix: 

Poo = — r^poo + r^pii , 
Pn = r^poo - Tjpn , 

poi = -iB\\p i - ^rpoi , (14) 

where the excitation rate Tf and relaxation rate are related to T\ and the 
equilibrium value (a z } by 1/Ti = T l + T T and (a z ) = T 1 (T i - T T ). 
A series of papers were devoted to the derivation and generalization of the 
Bloch equations (Wangsness and Bloch, 1953; Bloch, 1957; Redfield, 1957). 
Below we will illustrate the derivation from the model (8) in several lim- 
its of the qubit's dynamics. We can anticipate two different regimes: The 
Hamiltonian-dominated regime where the dissipative effects are slow com- 
pared to the Larmor precession. In this case it is convenient to describe the 
dynamics in the eigenbasis of the spin's Hamiltonian. The other, dissipation- 
dominated regime, arises, when the total magnetic field is weak. Then the 
dissipation dominates, and the preferred eigenbasis is that of the dissipative 
part of the Hamiltonian (oc a z ). 

3.2. Golden Rule and the Bloch equations 

In the eigenbasis of the spin (qubit) the Hamiltonian (8) reads 
1 1 

H = --AE a z - -X (cost/ a z - sin?? a x ) + H hath , (15) 
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where AE = y/B 2 . + B 2 and tan 77 = B x /B z . We denote the ground and 
excited states of the free qubit by |0) and 1 1), respectively. The coupling to the 
bath is a sum of a transverse (oc sin 77) and a longitudinal (oc cos rj) term. Only 
the transverse part can cause spin flips. In the weak-noise limit we consider 
X as a perturbation and apply Fermi's golden rule to obtain the relaxation 
rate, 1^ = r^^o), and excitation rate, Tf = ]T|o) >.|i) . E.g., for we obtain 

Ft = 2 ^^p l j:p t \mX\f)\^(E t + AE-E f ) 



h ■ 

2ir sin 2 77 

~h 4 

sin 2 7] f,^„ /,, vwvl , AiAE 



E« (i\x\f)(f\x\i) ^nJ dteq(Et+AE ~ Ef) 



Ah 2 



J dtJ^Pi (i\X(t)X\i) 



sin2 V / V 2 v n ~ 

= -4^- \ u>=AE/ft) ■ ( 16 ) 

Here |i) and |/) are the initial and the final states of the bath and pi is the 
probability for the bath to be initially in the state Similarly, we obtain 

r T = ( x I=-ae/h) ■ ( 17 ) 

For the relaxation time we thus find 

k = Tl + Ft = ^ SX(UJ = AE/k) ' (18) 
and for the equilibrium magnetization 

r,-r T A x (uj = AE/h) , AE 

The golden rule can also be used to evaluate the dephasing time. Here we skip 
the derivation since it will be performed in Section 3.5 in the framework of 
the Bloch-Redfield approximation. The result is 

T 2 X = ±Tf 1 + ^ cog2 V Sx(co = 0) . (20) 

Note that the first term oc S(u = AE/h) in Eq. (20) is the result of the 
transverse noise, which involves transitions between the qubit's states with 
the energy change of AE. The second term oc S{lo = 0) is associated 
with the longitudinal noise, which does not flip the spin and therefore in- 
volves only transitions in the bath without energy transfer. It produces a 
random contribution to energy differences, and hence to the accumulated 
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P s (0= P S (0) 



• = cosr| a z - sinr) a x 

Figure 2. Diagrammatic representation of the propagator IL 

phase difference. This contribution to the dephasing rate is sometimes called 
the "pure-dephasing" rate, 1/7%, so that T^ 1 = \ Tf 1 + (2J) _1 . 

3.3. Diagrammatic technique 

We proceed with a systematic perturbative analysis of the time evolution of 
the qubit's density matrix based on the Keldysh diagrammatic technique. We 
derive the Bloch equations and the corresponding dissipative times within 
the so-called Bloch-Redfield approximation, which may be thought of as a 
generalized golden rule. The systematic extensions allow us to analyze effects 
of higher orders in the noise and of long-time noise correlations. We will also 
determine the limits of applicability of the Bloch equations. 
First we briefly sketch the derivation of the master equation for the evolution 
of the density matrix, based on the formalism developed by Schoeller and 
Schon (1994). With the interaction Hamiltonian 

H int = ~\ x ( cos V °z ~ sin rj a x ) (21) 

the time evolution of the density matrix in the interaction representation fol- 
lows from 

t t 

-If H int (t')dt> | / H int (t')dt' 

p(t) = Te *o p(t ) fe <<> . (22) 

Assuming that initially, at t = to, the density matrix is factorized, p(to) = 
Ps(to) <8> Pbath, we trace Eq. (22) over the bath degrees of freedom and obtain 
a perturbative expansion for the propagator II of the spin's density matrix 
p s (t) = H(t,to)p s (to). On one hand, the assumption of a factorized initial 
condition allows us to introduce the characteristics of the bath, e.g. its tem- 
perature. On the other hand, it may introduce artifacts. In general we find 
for p s a non-Markovian (nonlocal in time) evolution equation such that the 
behavior at a time t involves knowledge of p s at earlier times. However, often 
the kernel of this integral equation has a limited range r c and does not depend 
on the initial spin state. Then the simple assumption about the initial density 
matrix should not influence the long-time evolution. 
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t' t t' 



= <X(t)X(t')> e -iH s (t'-t") e -iH s (t-f 

r \ / 



-• • •- 



L- t" f t 

= <X(t')X(t)> t > t' > t' 



t' t t' 

Figure 3. Basic elements of the diagrammatic technique. 



By tracing out the bath degrees of freedom in Eq. (22) and expanding in H int , 
one obtains a perturbative series. A typical contribution to II is shown in 
Fig. 2. The upper line contains vertices (the terms H mt ) from the first, time- 
ordered exponent in Eq. (22), and the lower one contains vertices from the 
other, inversely time-ordered exponent. As one can see from Eq. (21), each 
vertex in the diagrams contains one bath operator X. The emerging averages 
of their products reduce, due to Wick's theorem, in a standard way to pair- 
wise averages, which are shown by the thin lines in the figure. In Fig. 3a the 
diagrammatic rules for the thin lines are sketched. The horizontal solid lines 
(see Fig. 3b) describe explicitly the evolution of the spin degree of freedom, 
each line corresponding to (a 2 x 2 operator) exp[±iH s (t — t')] with different 
signs for the upper and lower lines. 

The propagator II satisfies the Dyson equation shown in Fig. 4. Taking the 
time derivative one arrives at the master (kinetic) equation 

t 

j t Ps{t) = l -[ Ps {t),H s ] + Jdt'X(t-t')p s (t') 

o 

t 

= L sPs {t) + J dt 1 £(* - t') Ps (t') . (23) 
o 

We introduced the Liouville operator L s p s = (i/h) [p s , H s ]. The self-energy 
£ is defined as the sum of all irreducible diagrams, i.e., diagrams that can not 
be divided by cutting the upper and lower solid lines at the same time (see 
an example in Fig. 5). It is an operator in the space of density matrices, i.e., 
a tensor of rank four. We denote its components by £ nn /.,_ mm / > implying the 
following action on density matrices: (Sp s ) nn / = ^ S nn / < _ w {p s ) mm >. 

mm 1 

3.4. Bloch-Redfield and the rotating-wave approximations 

Under certain conditions Eq. (23) may be approximated by a simpler equation 
which is local in time. Indeed, the self-energy (which in lowest order is pro- 
portional to the noise correlator) in the last term on the right hand side decays 
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n 



n 



(0) 



n 



n 



(0) 



Figure 4. The Dyson equation for the propagator II. 




t' t 
Figure 5. A (second-order) contribution to E nn /^ mm / (t — t'). 

on time differences t — t' of the order of r c . If this time is much shorter than 
the dissipative dephasing and relaxation times we may replace the density 
matrix in the integral from t' to t by the leading term 



Ps(t') 



-L B (t-f) 



Ps(t) , 



and obtain a Markovian evolution equation 



d 



—p s (t) = L sPs (t) + T Ps (t). 



(24) 



(25) 



It involves the 'Bloch-Redfield tensor' 



f= J dtt(t) 



,-L a t 



(26) 



To simplify the Fourier analysis here and in the following, the definition of 
the self-energy is extended to negative times with S(i < 0) = 0. In the 
eigenbasis of H s the Liouvillian L s is diagonal, [L s ] mm i^ mm / = i(E m — 
E m i)h, and the components of the Bloch-Redfield tensor are related to the 



self-energy: V 



{fttO = E m — E m i). 



To verify the validity of the Bloch-Redfield approximation one should check 
whether the dissipative times T\, T2 in the resulting Bloch equations are in- 
deed much longer than the correlation time r c of the integrand in Eq. (26) 
[more precisely, whether the integral is dominated by times shorter than T\ 
and T 2 ]. 2 



Note that due to differences in the oscillatory factors different components of the matrix 
integral in Eq. (26) can be dominated by different time scales. Within the RWA the influence 
of dissipative terms on the nn'-component of the density matrix is dominated by the behavior 
of the bath spectral power close to hu> = E n — E n i . One finds that the validity of the Bloch- 
Redfield approximation requires that the spectrum does not change much on the scale given 
by the respective component of f . 
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In certain cases the equations of motion can be simplified further by em- 
ploying the rotating-wave approximation (RWA). It is based on the following 
considerations: if the dissipation is weak and the dynamics of the components 
of the density matrix close to the unperturbed dynamics, the spectral weight 
of Pnn'(t) is located in the vicinity of to = (E n — E n >)/h = u^„/. The rhs 
of the nn'-component of the equation of motion (25) contains contributions 
from all components mm' of the density matrix. However, if w^ m ; / w° n / 
this contribution averages out fast (on the time scale ft/(u;^ m / — and 
can be neglected. I.e., one truncates the terms T nn r^ mm rp mm i on the rhs of 
the equation of motion: for the diagonal elements p nn one keeps only the 
diagonal entries, while for the off-diagonal elements, e.g. poi> one keeps only 
the contribution from the same entry, e.g. the component Toi^oi- The latter 
are, in general, complex numbers, with a real part describing dephasing and 
an imaginary part corresponding to a renormalization of the energy. 
If the system is subject to a strong pulse, Eq. (25) may not be sufficient. 
Typically the initial evolution for a certain period of time ~ r c needs to be de- 
scribed by other methods. However, the following time evolution is governed 
by (25). We will discuss examples in Section 5.3. Another situation, where 
Eq. (25) may be not sufficient is the case where the noise shows long-range 
correlations (cf. the discussion of the 1/f noise in Section 5.1). 



3.5. Ti AND T 2 

As an example and demonstration of the diagrammatic expansion we rederive 
the expressions for T\ and Ti, obtained before by golden rule arguments in 
Section 3.2. The lowest-order contributions to the rate roo<-n = ~£\ are 
shown in Fig. 6. We obtain 



sin 2 rj 



{X(t)X(t'))er, AE ^ + c.c] , (27) 



and thus, from Eq. (26) and the fact that (£ s )ii<-ii = 0' we § et 



r. 



Re 



00<-ll 



i sin r] 



S o^ii(^ = 0) = 
dJ {Xl,) 



2vr lo' - (AE/h) - iO 
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Figure 7. Leading-order contributions to rni<-oi . 



—^T ( X u=AE/h) ■ ( 28 ) 

In the same way one can obtain the excitation rate Tn^oo = If, and hence 
Ti, as well as the rates r o-oo = -Tn^oo and Tn^n = -Too-ii- 
As for the dephasing time T2, the relevant lowest-order contributions to Toi^oi 
are shown in Fig. 7. The upper row of the diagrams shows the contribution 
of the transverse fluctuations Eq^^ and the lower row shows that of the 

longitudinal fluctuations Sqi- oi' so that Soi^oi = ^oi— 01 + ^01—01 • 
The transverse and longitudinal contributions give: 

^ 01 (t-t') = - S -^S(t-t>), (29) 

4^ 01 (t - t') = - C ^p- S(t - t>) ei AE ^ . (30) 

Noting that (£ s )oi— 01 = i^-E/h we obtain from Eq. (26) 

, _ i sin 2 rj f duj S x (uj) 
1 01-01 " ~^-J 2^ u + (AE/h)-iO ' (31) 

II _ icos 2 r] f da; Sxjuj) 

1 01-01 - -p-y 5 



and Toi^oi = ^01— 01 +^01— 01 • The rea ^ °f these entries of the Redfield 
tensor, ReToi^oi = — > reproduces the dephasing rate (20). 

3.6. Energy renormalization 

The analysis of the preceding subsection also shows that the qubit's level 
splitting is renormalized by the bath. This effect is the analogue of the Lamb 
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shift, i.e. the renormalization of the level splitting in atoms due to the coupling 
to the electromagnetic radiation. As can be seen from Eq. (25) an imaginary 
part of Toi^oi gives rise to such a renormalization. Only the transverse noise 
contributes to this renormalization. From Eq. (31) we obtain 

Imro^oi = --20T- P J ^ - _ {AE/h) ■ (33) 

As an example we consider an oscillator bath with Ohmic spectral density 
up to a high-frequency cut-off, J(uj) = 2iraftwQ(uj c — u), at low tempera- 
ture k B T <C AE where S(AE) = J(AE) coth « J(AE), which is 
coupled purely transversely, rj = tt/2. Then the integral in (33) is dominated 
by high frequencies Huj > AE and gives the following correction to the 
energy splitting: AE — ► AE(1 — a ln(Huj c / AE)). This result indicates that 
the small parameter of the perturbative expansion is a ln(Hu c / AE), rather 
than a. Also for general bath spectra the high-frequency part of the spectrum 
typically leads to a renormalized Hamiltonian. The effect can be calculated 
perturbatively as shown above or in the frame of a renormalization-group 
approach, which accounts for all orders in the small parameter. It turns out 
that in such a situation only the component of the field perpendicular to the 
direction of the fluctuating field is affected by the renormalization (Leggett 
et al., 1987; Weiss, 1999). 

The effect of low-frequency fluctuations can be treated adiabatically : They do 
not induce transitions between the qubit states, but they change their energy 
splitting since the fluctuations of the transverse component of the 'magnetic 
field' X 1 - increase the average magnitude of the field. 
Even if the bare Hamiltonian is time-dependent the contribution of the high- 
frequency modes {ftw S> AE) of the environment can be absorbed in a 
renormalization of the Hamiltonian. In contrast, the analysis of the effect of 
low-frequency modes depends on further details. 



4. Applications 

4.1. Ohmic dissipation in real devices 

For an Ohmic environment we obtain from Eqs. (18) and (20) 

! , AE AE 

T{ 1 = ttq sin 2 v — coth ^-^ , (34) 

T 2 l = - Tf 1 + vra cos 2 r, ^ . (35) 

_■ lb 

For the charge qubit the coupling constant a is given by Eq. (12). The typ- 
ical impedance of the control line is R rs 50 £1. Since Rq « 6.5 kJ7 we 
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obtain a « 10™ 2 (Ct/Cj) 2 . With the typical values of the capacitances Cj « 
10~ 16 F and C g w 10~ 18 F their ratio can be made as small as 1(T 2 . Thus a 
can be as small as 10~ 6 . At temperature T ps 10 mK and for AE/k-Q ~ IK, 
we estimate the pure dephasing time T| ~ 100 /xs (for 77 = 0) and the 
relaxation time Ti « 1 fis (for 77 = 7r/2). In recent experiments (Vion 
et al., 2002; Chiorescu et al., 2003) similar values of T\ have been observed, 
however T2 was much shorter. This has been attributed to the contribution of 
the 1// noise (i.e. dominated by low-frequencies), which we discuss below. 

4.2. Validity of the Bloch-Redfield approximation in the 
Ohmic case 

As mentioned above, after having solved the Bloch equations one should 
verify the validity of the Bloch-Redfield approximation by checking the self- 
consistency. Let us consider the example of a weakly coupled Ohmic envi- 
ronment, a <C 1. The Bloch-Redfield approximation is justified provided the 
integral in Eq. (26) converges on a time scale short compared to the appropri- 
ate dissipative time T\ or T2. 

(i) For the evaluation of Tf 1 oc S(AE/h) sin 2 77 consistency requires that the 
noise power in the vicinity of uj = AE/h varies smoothly on the scale T-j -1 . 
For the Ohmic bath this condition is always satisfied: in the limit AE 3> ksT 
it reduces to aAE <C AE and in the opposite limit to aT <C T. 

(ii) For the evaluation of T^ 1 = ^Tf 1 + (T 2 *) -1 oc %S(AE) sin 2 r] + 
S(0) cos 2 7] the condition may differ only if the rate is dominated by the 
second term, T|. Then the condition is that the noise power varies weakly 
on the scale of (TI)™ 1 ~ aT{kB,/K) "C T (Ice/ ft) around zero frequency, 
which is always true. 

We conclude that for a <C 1 the Bloch-Redfield approximation is valid. 



4.3. Environment-dominated dynamics 



Consider a Hamiltonian with purely transverse coupling, rj 
H 



--A£ a z + -X a x + #bath 



For this model the Bloch-Redfield tensor becomes: 



~dt 



Here r„ 



/ Poo ^ 


( 


Pll 




P01 




V Pw ) 


\ 



/-r t 



-r 








iAE-Y 






-iAE -r*J 



vr/2: 



( Poo \ 

pll 

P01 

V "10 J 



(36) 



(37) 



. If the two-level system is almost degenerate, AE <C 
KTtp, the RWA is not valid, and one can not neglect the off-diagonal entries 



01<-01 
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(obtained from diagrams similar to those in Fig. 6). Notice that the degener- 
acy condition, AE <C T^, may still be consistent with the validity condition 
of the Bloch-Redfield approximation, l/r c S> T^. For instance, for an Ohmic 
bath in the weak-coupling limit, a <C 1, one finds l/r c ~ (k B /h)T » 
a(k B /h)T ~ r v (the latter equality is valid provided AE < k B T). 
In this regime, AE <C hT^, Eq. (37) has a very slowly decaying eigenmode, 
with p oc a x , with the relaxation rate (AE) 2 /2T V (T^ is almost real). This 
phenomenon is called the quantum Zeno effect (Misra and Sudarshan, 1977). 
The environment "observes" a x and prevents it from relaxing. The effect has 
the remarkable property that the stronger the dephasing, Y^, the slower the 
relaxation. 

In general, in the dissipation-dominated regime it is more convenient to treat 
the spin Hamiltonian as a perturbation on top of the dissipative dynamics. 

4.4. Time-dependent driving: Rabi oscillations 

For a general time-dependent magnetic field B(t) the description of the dissi- 
pative dynamics is more complicated. Here we discuss an example for which 
the Bloch-Redfield approximation does not apply in the laboratory frame 
but can be used in the rotating frame (or, equivalently, in the interaction 
representation). 

Consider a situation with the purely longitudinal coupling, 

H = ~AE a z - -X a z + H hath , (38) 

for which Eqs. (18) and (20) give Tf 1 = and T 2 l = (l/2h 2 ) S x (lo = 0). 
Let us now apply a resonant transverse driving field: 

H = ~7^AE a z - (cos u)ta x - smuta y ) - a z + H hath . (39) 

The subscript R indicates that this field can induce the so-called Rabi os- 
cillations between the levels. As follows from the discussion in Section 3.4, 
one can still use the Bloch equations as long as the driving is weak enough, 
^RTc *C 1. For a stronger driving these equations may fail. However, the 
transformation to the rotating frame [H — > H = hUU + UHU^, with U = 
exp (—iu>a z t/2) and Huj = AE] makes the Hamiltonian time-independent: 

H = a x -^Xa z + H hath . (40) 

The first term on the rhs of Eq. (40) induces coherent oscillations in the 
rotating frame, with the spin performing an oscillatory motion between the 
up- and down-states. Furthermore, using Eqs. (18) and (20), we find that the 
relaxation of the x-component of the spin in the rotating frame occurs on the 
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time scale T-f 1 = (\/2h 2 ) Sx{oj = Or), and the decay of the transverse (y- 
and z-) components on the time scale T^ 1 = l/(2Ti). One may verify that in 
the weak-driving limit (Ort c <C 1) this result in consistent with the standard 
Bloch equations in the laboratory frame. 

We summarize, strong or fast variations of the Hamiltonian (e.g., the field 
B(t)) may invalidate the Bloch equations, but in certain situations one may 
still reduce the problem to another one for which these equation can be ap- 
plied. The resulting dissipative times T\ , T2 may, however, differ from those 
in the original Bloch equations. 



5. Beyond the Bloch-Redfield (golden-rule) approximation 

5.1. 1/F NOISE 

Several experiments with Josephson circuits have revealed at low frequencies 
the presence of l/f noise. While the origin of this noise may be differ- 
ent in different circuits, it appears that in several cases (especially in the 
charge devices) it derives from "background charge fluctuations". This noise 
is usually presented as an effective noise of the gate charge (see Eq.(5)), i.e., 
Sq s (uj) = a.i/fe 2 /\ui\. Recent experiments (cf. Nakamura et al. (2002) and 
references therein) yield at relevant temperatures a^jf ~ 10~ 7 -10 -6 . We 
translate this noise into fluctuations of the effective magnetic field X with 

E 2 

S x (u) = -±tf , (41) 

\LU\ 

and by comparison with Eq. (9) find Eyj = 4Ec^/aT/J- 
The Bloch-Redfield result (20) predicts a divergent rate T 2 _1 — > 00, unless 
77 = 7r/2, and clearly fails in this case. This can be attributed to the very slow 
(logarithmic) time decay of the function Sx(t). It is easy, though, to describe 
the dephasing due to 1// noise rigorously, if it couples purely longitudinally 
(77 = 0) and is Gaussian distributed. Treating X as a classical variable one 
obtains (Cottet et al., 2001) 

/ -| / dt'X(t')\ -T^fdt' J dt"(X{t')X{t")) 
(cr + (t)) oc ( e ) = e » 



dco „ . , sin 2 (u;£/2) 

— — ij w 

2vr XK ' ( u /2y 



exp(-7^ I^S x ( U )^p-\ . (42) 



This result coincides with the real part of Eq. (55) obtained in a quantum 
description. For a regular low-frequency spectrum one may replace the last 
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fraction in (42) by 2tt5(lo) t and obtain the golden-rule result for the dephas- 
ing time. For the 1/f noise one obtains with logarithmic accuracy (Cottet 
et al., 2001) 

(E 2 \ 
-^t 2 |ln^ ir |j , (43) 

where a; ir is the infrared cutoff frequency for the 1/f noise. 3 From this decay 
law one can deduce the dephasing time scale, 

(44) 

Using experimental parameters for an estimate one obtains T 2 in the range 
from fractions to a few nanoseconds, which is close to the dephasing time 
observed in the experiments away from the special point rj = tt/2. We also 
see that this source of dephasing dominates over the Ohmic noise. 

5.2. Contribution of higher orders 

In this section we consider a further situation where the lowest-order Bloch 
equations are not sufficient. The dissipative times (18, 20) in the Bloch equa- 
tions are dominated by the noise power at specific frequencies (0 and AE/h). 
As an example, consider the purely transverse case, rj = tt/2, where the first- 
order Bloch-Redfield formulas give Tf 1 = 2T 2 X = ^S x {u = ^p). In 
higher orders of the perturbative expansion (cf. Section 3.3) the noise power 
in other frequency ranges may be involved. If the noise power in these ranges 
is much stronger, higher-order contributions may compete with the lowest- 
order. An important example is the case of 1// noise. Here the strong noise 
at low frequencies contributes in higher order to the dissipative rates. 

To illustrate the importance of higher orders, consider the contribution to the 

(2) 

second-order self-energy Sq/ < _ 01 shown in Fig. 8a, 

t Tl 



v 



(45) 

iAE 1 _ \ 

The factor e r ^ 1 T2) is due to the difference in energy of the state on the 
upper (|0)) and lower (|1)) lines in the time interval between n to t 2 . The 
Fourier transform of Eq. (45) reads 



AV( 2 ) f,,\ ~ l f duJ l duJ 2 <v2u Y 2> 



3 In an experiment with averaging over repetitive measurements the infrared cutoff may be 
determined by the time interval t av over which the averaging is performed (Vion et al., 2002). 
In a spin-echo experiment it may be set by the echo frequency (Nakamura et al., 2002). 
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t 
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t 



Figure 8. Second-order contributions. One can obtain all 16 second-order diagrams by 
shifting an even number of the vertices in a or b to the lower line. 



x \ (46) 

(uj - ui + i0) 2 {uj -U1-W2 + ^E/h + »0) 

Comparison of Eq. (46) and Fig. 8a elucidates the diagrammatic rules: the 
bath lines carry the frequencies uj\ and uj2, while pairs of horizontal solid 
lines contribute factors of the type 1/ (a; — X) a> n — iL B ), where the summa- 
tion is over all the bath lines in the interval and the matrix element of L s is 
determined by the spin states on the upper and lower lines. 
There are in total 16 second-order diagrams (see Fig. 8), which contribute to 
Sqi<-oi- They add up to 



,(2) 

J 01^01 



00 

/ 



-2i 



{2hY J (2tt) 5 



■ + 



x 



—00 
1 



uj — cji + iO uj — uj2 + iO 

\xl 1 ){xl 2 ) + {xl^){x^ 2 ) 

uj — UJ\ 



uj — uj\ + iO uj — uj\ — UJ2 + AE/h + iO 

, {xl 1 ){xl u)2 ) + {x^ 1 ){xl 2 ) 

uj — uj\ — UJ2 — AE/h + iO 



(47) 



We evaluate the behavior of the self-energy in the vicinity of the level split- 
ting, EqiLoiy^ = -AE/h - uj' + iO), where uj' < AE/h. If the integral 
(47) is dominated by low frequencies, we obtain: 



AE 



-uj' + io) 

00 

1 f dv 



8h 2 AE 2 J 2ir 

—00 



= - W S X2 (uj') 



(Xl +u)l )(X 2 v ) + (Xl_ u ,)(X 2 v ) 



(48) 



where Sx 2 is the symmetrized correlator of X2 = X 2 / (2AE). 
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To understand the origin of this result (48) consider the effect of the low- 
frequency fluctuations (which, we assumed, dominate the dephasing). Start- 
ing from the Hamiltonian (15) with the transverse coupling to the noise source, 
r) = it/2, we see that in the adiabatic limit of very slow fluctuations the 
dynamics reduces to an accumulation of the relative phase between the eigen- 
states, JdtAE x (t)/h, where AE(X) = ^ AE 2 + X 2 « AE+X 2 /(2AE). 
Thus the problem reduces to the longitudinal case. The calculation above 
confirms this reduction in the lowest order. More precisely, the second-order 

contribution to Sqi^-oi ^ or tne transversely coupled X coincides with the 
lowest-order contribution for the longitudinally coupled X2. 
Notice also that the fluctuating field X2 has a nonzero average which provides 
a (low-frequency contribution to the) renormalization of the level splitting 
(cf. Section 3.6). 

Generally speaking, the result (48) does not allow one to evaluate the de- 
phasing decay law but may be used for qualitative estimates. Consider the 
example of the 1/f noise. In this case the low-frequency behavior of the 
second-order self-energy is singular and it may compete with the first-order 
contribution. Higher-order contributions require further analysis (both in the 
linear transverse and quadratic longitudinal cases). However, one can show 
that at high enough frequencies the higher-order contributions to £ may be 
neglected. Thus, the result (48) permits an analysis of the decay law at short 
times. Extrapolating this law to longer times, one can estimate the typical 
dephasing time. Qualitatively, this amounts to solving the equation T^ 1 ~ 

— £ 01 ;_ 01 (u; = T 2 ). For 1/f spectrum (41) we obtain 

E-\lf 1 \uj\ 

SxM-j^Ly-M—, (49) 
(AE) Z \uj\ uj- n 

i.e., this noise power is again of the 1/f type. This gives the estimate T 2 _1 ~ 
{El /f /hAE)\n(E 2 l/f /hAEio ir ). 

5.3. Effect of strong pulses 

In this section we discuss the effects of a strong pulse of the magnetic field 
on the dissipative dynamics of the spin. The derivation of the Bloch equa- 
tions (cf. the discussion in Section 3.4) relied on the slowness of the spin's 
dynamics, which may be violated if a strong or a fast-oscillating field is 
applied. In this case for a period of order r c a different description beyond 
the Bloch-Redfield approximation is needed. 

As an example we consider the 'initial-state preparation' of a gedanken ex- 
periment designed to observe the effects of dephasing on the spin precession. 
Specifically, we assume that initially the spin was kept at low temperatures 
in a field B z z and has relaxed to the ground state | j). Then it is rotated by a 
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strong short 7r/2-pulse of the transverse field and starts a free precession in 
the field B z z. Under the influence of the bath the relative phase between the 
spin-up and -down states will not only have a contribution due to the Zeeman 
splitting, which grows linearly in time, but also a random contribution from 
the noise, which destroys the phase coherence. For the analysis below we 
assume the Ohmic spectral density. 

The effect of the pulse on the spin's dynamics, as compared to the Bloch 
equations, is an additional unitary evolution and dephasing. The former effect 
is an extra phase acquired by the spin due to the initial polarization of the 
bath. The latter is the dephasing due to the initial partial disentanglement of 
the bath and the spin (which amounts to a partial factorization of the initial 
density matrix), as we illustrate below. 

During the relaxation of the spin to the ground state ||) the oscillators have 
also relaxed to the corresponding ground state. I.e., the initial state of the 
whole system is given by the density matrix po = |T) (t| <8> P\ with oc 
exp(— H^/k&T) and = i^bath — X/2. In this state the bath operator 
X provides a finite average field. Its expectation value (X) can be calcu- 
lated as a linear response of the bath to the perturbation H VCTt = — fX 
with / = l/2, i.e., {X) = x (w = 0)/ = = 0)/2. For an Ohmic 

spectrum, x"( w ) = ^Trahuj up to a cutoff uj c , we have (X) = x'(oj = 
0)/2 = f(dv ~x" {v)) / (2itv) ~ ahuj c . We note that this expectation value 
depends strongly on the high-frequency cutoff lo c of the bath spectrum. The 
field produced by each bath mode averages out after its typical oscillation 
period; still during a short initial period of length ~ uj^ 1 it affects the spin 
dynamics. Depending on the spectral density of the bath, the resulting effect 
may be dominated by low or high frequencies. 

Now we apply a 7r/2-pulse to rotate the spin to a position perpendicular to the 
original state, i.e. into the superposition ( | T) + 1 1) ) / V%- If the pulse duration is 
very short, the state of the bath does not change during the pulse. If the pulse 
takes a finite time, u;" 1 , the bath oscillators partially adjust to the changing 
spin state. For the high-frequency oscillators, uj 3> w p , the spin evolution 
during the pulse (and afterwards) is adiabatic; they just 'dress' the spin, renor- 
malizing its parameters, e.g., its g-factor (Shnirman et al., 2002). In contrast, 
the slow oscillators, with lo <C lo p , do not have time to react on this short time 
interval and start to change their states only after the pulse. They react dif- 
ferently for the spin-up and spin-down components, thus contributing to the 
dephasing. Neglecting, for a qualitative discussion, the oscillators at interme- 
diate frequencies, we find the system in the state po = \ T + 1^ + I ®f>\- 

Here the state f ^ denotes the state of the dressed spin, i.e. the product state of 
1 1 ) and the equilibrium states of high frequency oscillators under the influence 
of the spin up (i.e. with the Hamiltonian H^), and similarly for the state J). 
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In short, qualitatively one may think of the system as a dressed spin coupled 
to a bath with the upper cutoff ~ uj p . 

After this initial preparation, the spin starts a precession in the field B z z. To 
describe the time evolution in the presence of the bath we split the Hamilto- 
nian (8) as 

H = -\B z a z - X -X{o z - 1) + # T (50) 

and extract explicitly the average (X) with the intention to treat the part 
proportional to 6X as a perturbation. Hence 

H = -\b z o z - \{X){a z - 1) - \bX{a z - 1) + F T . (51) 

At t = the spin is subject to the additional field (X). We expect that 
this field relaxes, on the time scale uj^ 1 for oscillators with frequency u. 
Nevertheless, before the field vanishes, it contributes to the phase evolution. 
Clearly, the spin-up and spin-down states are decoupled, which allows for an 
exact solution. One may use the so called polaron transformation (cf. Leggett 
et al. (1987), Weiss (1999)), or calculate the time evolution of the off-diagonal 
element of the qubit's density matrix directly: 

<<7+(t)> = e-* B ** Tr [a+ S(t,0) Po S\t,0)] , (52) 

where 

i 

^JSX(t')(a z -l)dt' 

S(t, 0) = Te o . (53) 

Obviously, the evolution operators S and do not flip the spin, while the 
operator a+ in Eq. (52) imposes a selection rule such that the spin is in the 
state ||) on the forward Keldysh contour (i.e., in S(t,0)) and in the state 
||) on the backward contour (i.e., in S^(t, 0)). Since (|| a z — 1 1|) = the 
contribution of the backward contour vanishes. From (|| a z — 1 1|) = —2 we 
obtain (o+(t)) = P(t)exp(-iB z t/h)(a + (0)), where 

P(t) = (Texp ^-^ J 5X(t')dt^j ) . (54) 

Due to the harmonic nature of the bath the equilibrium correlation functions 
of 5X for a bath with the Hamiltonian coincide with those for X and 
#bath- Since the fluctuations of 5X are Gaussian, we find 

t t 



In P(t) = J J dt 1 dt 2 G c (t 1 -t 2 ) 



o o 



i f du S m 2 (u;t/2) 
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In equilibrium the Green function G c (t - 1') = -i(T8X(t)6X(t')) is related 
to the noise by the fluctuation-dissipation theorem: 

G C M = -h (x'H + ix» ™th |0 = -h x 'H - iS(u) . (56) 

If we substituted sin 2 (cot/ '2) /(uj/2) 2 -> 27r<5(u;)i in Eq. (55) we would 
obtain 

lnP(t) = ~w Sx(uj = 0)t+ ih x ' iuj = 0)t - (57) 

The first term in Eq. (57) yields the standard longitudinal dephasing rate 
1/T| = 2irak-QT /h, while the second term exactly cancels the effect of the 
additional magnetic field (X). Thus, this approximation predicts a vanishing 
dephasing rate at T = 0, and there is no net effect of the initial preparation. 
However, the substitution above is justified only if G c (uj) has no structure at 
low frequencies oj ~ 1/1 This condition is violated, for an Ohmic bath, for 
not too long times, still small compared to the inverse temperature uj p 1 < 
t < h/k-gT (cf. e.g., Shnirman et al. (2002)). In this time range we find 
Re [lnP(i)] — 2a\ ln(cj p t)|, which implies a power-law decay 

\(a + (t))\ = (u p tr 2a \(a + (0))\ . (58) 

Thus even at T = 0, when Sx(oJ = 0) = 0, the off-diagonal elements of 
the density matrix decay in time. The reason is the excitation of the low- 
frequency (lo < Up) modes of the environment by the initial strong pulse. 
Only the excited modes contribute to the dephasing. From this viewpoint, the 
effect of the initial pulse is similar to heating. 

Finally, the extra phase contribution 5$>(t) = Im [In P(t)] — (X) t/h is 
n J 2ir uj z 

In deriving Eq. (59) we have used the analyticity of \ m tne upper half- 
plane. For the Ohmic bath, J$ — ► — ira for t S> u p 1 . Thus the extra phase 
is acquired during a short period of time of order uj~ 1 . One may attribute 
this extra phase to the additional magnetic field (X) produced by the excited 
modes in the bath. Note also that the time, over which the extra phase is 
acquired, may depend on the power spectrum of the noise. 
After a period of free precession a measurement of the final spin state is to 
be performed in order to monitor the effect of dissipation. A 'low-energy 
measurement', which does not break the adiabaticity condition for the high- 
frequency modes, reads out the state of the dressed spin. 
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5.4. 'Optimal' operation points and nonlinear coupling 

In the presence of 1/ f noise, its longitudinal component typically dominates 
dephasing (see Section 5.1). By adjusting the system control parameters to 
a point where the longitudinal coupling to the low-frequency noise is tuned 
to zero one should be able to increase the coherence time. This increase was 
observed, indeed, in recent experiments (Vion et al, 2002). The qubit in this 
experiment was subject to several electromagnetic noise sources, including 
charge and flux noise. When tuned to the optimal operation point the re- 
maining coupling to the noise is linear transverse for one part, X, of the 
noise sources and quadratic longitudinal for the rest, Y. In the mentioned 
experiments X refers to the charge noise, while Y describes flux noise. The 
resulting Hamiltonian reads 

H = ~(AEa z + Xa x + Y 2 a z ) . (60) 

Here, similar to the discussion after Eq. (48), we introduced Y 2 = Y 2 / (2Ey). 
The energy scale Ey characterizes the strength of the quadratic coupling to 
the noise source Y . For the setup of Vion et al. (2002) Ey is of the order of 
the qubit's Josephson energy. 

For noise spectra with a low-frequency singularity (e.g., 1 / /), both the trans- 
verse linear and the longitudinal quadratic coupling can induce strong de- 
phasing. As discussed in Section 5.2, for a transverse coupling to a source of 
low-frequency noise it may be necessary to account for higher-order contri- 
butions in the calculation of the dissipative rates. In particular, we have seen 
that, at least in the second order, the contribution of the transverse noise X 
reduces to that of the longitudinally coupled field X 2 = X 2 /(2AE). Thus, 
both the linear transverse coupling and the quadratic longitudinal coupling in 
(60) have similar effects on the dephasing. In this approximation, the effective 
low-frequency Hamiltonian reads 

H = ~AEa z -^(X 2 + Y 2 )a z . (61) 

We neglected the correlations between X and Y, and can estimate their ef- 
fects independently: (a + (t)) oc Px 2 {t)Py 2 (t), where the functions Px 2 and 
Py 2 characterize the dephasing due to the respective noise sources. As an 
example we estimate Px 2 (t). The noise power of the variable X 2 is given 
in Eq. (49). Treating X 2 as Gaussian distributed we obtain by substituting 
Eq. (49) in Eq. (42) the decay law (see Table I) 

(E^ \ 
-^^i 2 m 2 (^)] (62) 

and an estimate for the dephasing rate, l/T 2 ~ (E 2 ^/hAE) ln[E 2 ^/ AELO ir ] 
Note, however, that the result (62) is valid only at short times t C T 2 *. 
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Table I. The decay law for the off-diagonal (in the qubit's eigenbasis) entries of the qubit's den- 
sity matrix. The Bloch-Redfield approximation predicts an exponential decay with dephasing time 
oc S' _1 (a; = 0). Three examples are listed where this time vanishes or diverges and a more careful 
analysis is required. The sign || indicates the longitudinal coupling (rj = 0). 





Bloch-Redfield 
Sections 3, 4 


V/. II 

S(oj) oc E 2 1/f /uj 
Section 5.1 


x\ iff, II 

Sx 2 (u) oc In ^-/u) 
Section 5.4 


Ohmic, || 
S(ui) oc acu 

Section 5.3 


poi oc 


e -t/T 2 


e -El /f t 2 \\*(u lt t)\ 


-!$U-t 2 ln 2 (^ ir i) 


M)- 2 «fort<^ 



At longer times t ~ T| the non-Gaussian nature of the fluctuations of X 2 
becomes important (Makhlin and Shnirman, 2003). These, however, are not 
expected to substantially modify the estimate for the dephasing time scale 
obtained from an extrapolation of (62). In conclusion, we note that for AE » 
Ei /f the dephasing due to 1/ f noise at the optimal point is substantially 
reduced as compared to the general case where X is coupled linearly to a z . 



6. Summary 

In this tutorial we have discussed the issue of dissipation in solid-state qubits 
due to the coupling to the environment. First, the electromagnetic environ- 
ment of the charge qubit was described and the coupling Hamiltonian was 
derived, and the relevant Caldeira-Leggett model was introduced. For weak 
dissipation the description by the Bloch equations is sufficient in many sit- 
uations. We have presented the calculation of the relaxation and dephasing 
times, T\ and T 2 , in the Bloch-Redfield approximation and analyzed their 
validity range. We have also considered examples where the Bloch-equation 
description is not sufficient. This includes the effects of the strong pulses and 
dephasing due to 1/ f noise. Some of the examples analyzed are illustrated in 
Table I. It includes situations in which the Bloch-Redfield approximation is 
insufficient since it would produce vanishing or diverging dephasing rates. 

The work is part of the EU 1ST Project SQUBIT and of the CFN (Center for 
Functional Nanostructures) which is supported by the DFG (German Science 
Foundation). Y.M. was supported by the Humboldt foundation, the BMBF, 
and the ZIP Programme of the German government. 
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Z(CO) 



Figure 9. The qubit's environment modeled by a transmission line. 

7. Appendix: derivation of the Hamiltonian 

In this section we present a way of introducing the environment following 
Ingold and Nazarov (1992), in which the bath is modeled by a transmission 
line. Even if a real device is not coupled to an LC-line or a resistor, in generic 
situations the effect of the fluctuations on the system (qubit) can be modeled 
by an LC-line. The transmission line is also equivalent to an oscillator bath 
model introduced by Caldeira and Leggett. The purpose of the following 
is to present a particular model which might clarify some ideas behind the 
Caldeira-Leggett approach. A similar analysis for a slightly different model 
was presented by Paladino et al. (2003). 

Consider a transmission line, or LC-line, shown in Fig. 9. It amounts to an 
impedance Z{w) 

z-\u) = y (-iuL a - -L- + oV 1 = y — ^ , 

U iooC a J ^ L a (uj 2 -u; 2 a + iO sign uj) 

(63) 

with uj a = 1/ \JL a C a (index a counts the LC elements of the transmission 
line). The infinitesimal resistance in each LC-element serves to regularize the 
model. One way to derive the Hamiltonian is to start by introducing the de- 
grees of freedom (the phase drops at circuit elements), writing the Lagrangian 
of the system (a combination of the charging and Josephson energies), and 
then performing the Legendre transformation. For convenience instead of the 
phase variables 6^ (k being any index) we use the dimensionful flux variables 
4>k = c &o^fe/2vr (defined as time integrals of respective voltage drops; here 
$0 = h/2e), to which we still refer as 'phases'. As independent variables 
one can use the phase on the Josephson junction 0j, the phase between the 
terminals of the impedance 4> an d the phases on the capacitors of the LC 
elements of the transmission line 4> a . Then the phases on the inductances of 
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the LC elements are given by <f> — 4> a , while the phase on the gate capacitor 
is given by the Kirchhof rule 4> g = —(4>j + </> + / V g dt'). The Lagrangian of 
the system reads 



L = 



where 



+ 



C g (fa + ^ + V g ) 2 



- Ej cos 



2^)+L z (<f>), 



(64) 



C„ 



2L n 



(65) 



From here it is straightforward to derive the Hamiltonian of the system, 



H = 



(2en 



where 



2Cj 



2a 



2a 



+ 



(<f>a - <j>? 

2L a 



(66) 



(67) 



Here (n,9), (q,(f>), (q a ,4> a ) are the pairs of conjugate variables, and q = 
—C g 4>g is the charge on the gate capacitor. We assumed that the charge of 
the island 2en is completely located on the plates of the capacitors Cj and 
C g , and neglected the so-called offset charge induced by parasitic voltages 
due, e.g., to extra charges in the substrate. Introducing the new charge q = 
q — C t V g we obtain 



(2en - C t V g f 
H= 20 



2e 

- Ej cos 9 - —nq + C t V g 
w 



> + 



^- + J ff z (0), (68) 



thus recovering Eq. (7). 

Indeed, the degrees of freedom q a , 4> a and q, 4> play the part of the oscilla- 
tor bath, its Hamiltonian being given by the last two terms in Eq. (68) (its 
eigenmodes are not given explicitly). Further, qj C t = 5V corresponds to the 
voltage fluctuations, and the first term in Eq. (68) coincides with the sum of 
the first and the last terms in Eq. (7). Notice that the last two terms in Eq. (68) 
describe a parallel connection of the impedance and an effective capacitance 
Ct. Accordingly - one can check this also by an explicit calculation - the pro- 
duced noise is given by Eq. (6), i.e., the bath 'corresponds' to the impedance 
Z t rather than Z. 

In this model one can also observe the result of an abrupt change of the gate 
voltage. It is most easily seen from Eqs. (66, 67) that the gate voltage is 
coupled to the qubit's degree of freedom n only via the damped mode (q, (f>). 
As a result, for a pure resistor Z(uj) = R changes of V g influence the qubit 
only after a delay time, of order RCt- In the form (68) the term proportional 
to the time derivative of V g is responsible for this delay. 
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